Efficient frequency conversion through absorptive bands of the nonlinear crystal 
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Two simultaneous three wave mixing processes are analyzed, where an input frequency is converted to an 
output frequency via an intermediate stage. By employing simultaneous phase-matching and an adiabatic mod- 
ulation of the nonlinear coupling strengths, the intermediate frequency is kept dark throughout the interaction, 
while obtaining high conversion efficiency. This feat is accomplished in a manner analogous to population trans- 
fer in atomic stimulated Raman adiabatic passage (STIRAP). Applications include conversion between remote 
frequencies, e.g. mid-IR to visible, and study of electronic crystal properties in the UV absorption band. 

PACS numbers: 42.65.Ky, 42.65.Sf 



I. INTRODUCTION 

Quadratic nonlinear optical media are commonly used to 
perform optical frequency conversion via three wave mixing 
(TWM) processes [IJ. In this way, laser frequencies which 
are not available by direct laser action are generated. Com- 
binations of such processes were utilized extensively in order 
to reach frequencies either far removed from a laser source or 
close to it Jj, 01- This task is not trivial, since in order for a 
TWM process to be efficient a phase-matching condition has 
to be fulfilled. For a single TWM process, a well known solu- 
tion is quasi phase matching (QPM) |4], e.g. by periodically 
modulating the sign of the nonlinear coefficient. Schemes also 
exist for simultaneous phase-matching of multiple processes 
JUI^, at the expense of conversion efficiency. Furthermore, 
the medium is required to be transparent at all participating 
frequencies, including any frequency generated as an inter- 
mediate stage before further conversion. 

Recently, analogies were shown between the dynamics of 
simultaneous TWM processes and those induced in three level 
atoms by electromagnetic fields. In one publication, an anal- 
ogy is utilized to perform second harmonic generation and 
sum frequency generation (SFG) simultaneously, such that the 
third harmonic is generated efficiently [7]. In this case, the 
dynamics is nonlinear and significant power was generated 
at the intermediate frequency (second harmonic). In another 
work the analogy was applied to the case of purely linear dy- 
namics, where two simultaneous sum or difference generation 
processes (DFG), or a combination of them, were performed 
under the undepleted pump approximation |8]. With this 
method, called adiabatic elimination, only negligible power 
is ever present at the intermediate frequency throughout the 
entire conversion process, so transparency at the intermedi- 
ate frequency is not required. However, this method explicitly 
requires that each one of the two processes has a large phase- 
mismatch. Due to inherently large phase-mismatches, high 
conversion efficiency is difficult to achieve, and requires high 
pump intensities or very long interaction lengths. 

Here, we propose a method to perform two simultaneous 
TWM processes efficiently and without significant generation 



of an intermediate frequency. In contrast to adiabatic elim- 
ination 1 8], the two processes are phase-matched throughout 
the entire interaction length. High efficiency can therefore be 
maintained even in the presence of significant absorption of 
the intermediate frequency wave. As a result, one of the non- 
linearity enhancement mechanisms which is associated with 
absorption resonances can be exploited. In this case, the anal- 
ogous atomic scenario is the so-called stimulated Raman adi- 
abatic passage (STIRAP) j^ fioll . We emphasize that here the 
dynamics are linear, in complete correspondence to STIRAP. 
It is this linearity which enables the elimination of the inter- 
mediate frequency, in contrast to the previous work by Longhi 
Jt!], which mimicked a nonlinear version of STIRAP, and con- 
sequently significant intermediate power was present (^ 20% 
of the input power). 

The insensitivity to intermediate frequency absorption 
makes possible processes in which the intermediate frequency 
is in the ultraviolet (UV) absorption band of a nonlinear crys- 
tal. This method can thus be used to probe the electronic prop- 
erties of the crystal in this spectral region, a task which was 
previously made difficult by absorption. 

This paper is organized as follows. In section|II] the theoret- 
ical model is presented and adiabatic interaction is analyzed. 
In section Uni two study cases, which demonstrate the main 
features of our method, are simulated numerically assuming 
ideal conditions. Section|IV]repeats the calculation of section 
|III]under conditions that are attainable with current QPM tech- 
nology, and discusses its limitations. In section [Vl we discuss 
fuither improvements of the QPM -based method and another 
way of achieving nonlinear optics STIRAP, as well as appli- 
cations of this technique. 



II. THEORETICAL MODEL AND ANALYSIS 

A. Theoretical model and atomic analogy 

In this model, two simultaneous TWM processes are con- 
sidered, such that the frequency generated by one process is 
further combined with another frequency in the other proces s. 
For example, a SFG process 0)2 = Oi + (Op, can be accom- 



panied by a DFG process 0)3 = 0)2 — (Op-, = 0)i + 0),,, — 0)^ 
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where 0)i is the input frequency, 0)2 is the intermediate fre- 
quency, 0)3 is the output frequency and 0)p, and (Op^ are the 
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first and second pump frequencies, respectively. We make 
the undepleted pump approximation: the two pump waves are 
taken to be much more intense than the other waves, thus they 
are negHgibly affected by the interaction. Furthermore, all 
beams are assumed to be plane waves. If one uses Gaussian 
beams, this requirement can be translated to a Rayleigh range 
Zr that satisfies zr > L/2, where L is the nonlinear medium's 
length. 

The corresponding coupled wave dynamics equations for 
this scenario are 
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where \ ~ [Ai,A2,AiY' is the state vector, in which xj/j = 
Aj (z) is the complex envelope of the amplitude of the wave 
with frequency coj, and 
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the coupling matrix 
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are the effective coupling 
and Ap, are the complex 



coefficients between the fields. A 
envelopes of the amplitudes of the pumps fields, x^^' is the 
second-order nonlinear coefficient of the material and c is the 
velocity of light. A^i = ki± kp^ — k2 and Ak2 = k2± kp^ — k^ 
are the phase mismatches of the two nonlinear processes, 
where kj — nj {(Oj) COj /c is the wavenumber of the wave 
with frequency (Oj. The undepleted pumps approximation 
dApJdz = dApJdz = is included implicitly. These 
equations can describe four different cases: two SFG pro- 
cesses, two DFG processes, SFG followed by DFG and DFG 
followed by SFG (see fig. [T]i. The difference between these 
processes is manifested in the choice of either the top or 
bottom sign in )f,y and Akj. The top sign corresponds to SFG, 
while the bottom sign corresponds to DFG. 

These equations are isomorphic to the dynamics equations 
of a three level atom interacting with two EM fields |9, 10]: 
the optical interaction takes place over space (the z axis in 
our notation), while the atomic interaction takes place over 
time. The complex amphtude of each field Aj corresponds 
to the probability amplitude of population in each atomic 
state, often denoted Uj. Each effective coupling coefficient 
Kij is the counterpart of the dipole coupling strength between 
atomic levels / and j, usually described using the Rabi fre- 
quency (?) — dij ■ e{t) /Ti, where d,] is the dipole mo- 
ment between levels / and j, e (f) is the induced EM field, 
and h is the reduced Planck's constant. Finally, the phase 
mismatches, Aki and Ak2, correspond to the detunings with 
the two atomic transitions, Ai and A2, respectively. We note 
one minor difference: due to symmetry considerations, in the 
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Figure 1. The four cases described by Eq. [Tj (a) two SFG processes 
(b) two DFG processes (c) SFG followed by DFG (d) DFG followed 
by SFG. 



three level atom £2,^ = £2*-,- when the two coupling lasers have 
the same peak intensity, while in the nonhnear optics case 

Kij — {^cofkj / co^ki^ Kjj 7^ fCy,-, as noted above. 

As explained and demonstrated in previous works ifTTI . 
these equations, which assume quasi-monochromatic laser 
beams for the optical case, are applicable to pulses as short 
as 1 ps, owing to the mismatch between the group velocities 
of the different interacting wavelengths. Shorter pulses can 
be stretchered before the interaction and compressed to their 
transformed limited duration after the interaction [12,1 . 



B. Analysis of interaction under adiabatic variation 

A brief description of atomic STIRAP will be given here, 
while a more rigorous analysis will be performed for the non- 
linear optics variant. In atomic STIRAP |9, 10], atomic popu- 
lation is transferred from an initial state to a final state, where 
there is no dipole coupling between these two states. Rather, 
each of these two states is dipole-coupled to an intermediate 
state. Interestingly, STIRAP achieves this population transfer 
without ever having significant population in the intermediate 
state. This feat is accomplished by adiabatically varying the 
Rabi frequencies £2,y (f ) via control over the temporal shape of 
the induced EM field e (f ). Specifically, a pulse that couples 
the intermediate and final state is introduced into the atomic 
system, followed by another pulse which couples the inter- 
mediate state with the initial state, such that these two pulses 
have some temporal overlap. Only when this counter-intuitive 
pulse order is used is STIRAP achieved and the intermediate 
state remains unpopulated. This absence of population in the 
intermediate state renders the population transfer insensitive 
to radiative losses which may result from the short life time 
of the intermediate state, hence in this context it is termed a 
'dark state'. 

Note that single-photon resonance Ai = A2 = is not re- 
quired for STIRAP; two-photon resonance Ai + A2 = is suf- 
ficient. However, for the remainder of this section we will 
consider the case of perfect phase-matching Aki = Ak2 = 0, 
analogous to the atomic single-photon resonance. 
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An analysis of nonlinear optics STIRAP will now be given 
in a manner completely analogous to the standard analysis 
given in the atomic case jgl nOR . The three eigenvalues of 
the coupling matrix M are 



JfO = , K-± = ±VK'12K'21 + K'23K'32 (3) 

and their corresponding normahzed eigenvectors are 
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where we have defined fQ ~ 



H2 



Kij and — 



VK-12K-21 + K-23K-32- 

Next, let us define the angle : 



K'32/ 



(5) 



\g()) can now be written as 




(6) 



This eigenvector has an important property: it has no con- 
tribution from the intermediate state, i.e. (^o|A2) — 0. It is 
this property which will enable energy to be transferred from 
Ai to A3 without going through A2, as will now be explained. 

We assume that initially all of the optical power is in Ai, 
and additionally | K"i2| , | K'2i I < |k"23| , |k'32|. so w 0. In this 
situation, out of the three eigenstates, only \go) has a nonzero 
term, thus the system can be said to be in this eigenstate. 
During the interaction, the coupling coefficients change adi- 
abatically, such that at its end |k'i2| , |k'2i| ^ |k'23Mk'32|. so 
9 « n/2. According to the adiabatic theorem, if the system 
is at an eigenstate and is subject only to adiabatic changes, it 
will remain at the same eigenstate, in this case \go). However, 
when 9 « k/2, the physical meaning of {go) is that all of the 
optical power is in A3. We therefore conclude that the system 
will remain in the |^o) eigenstate throughout the interaction, 
going from the initial state [1,0,0] (all power in Ai) to the 
final state [0,0,-1]^ (all power in A3). Since \go) contains 
no component of the intermediate frequency amplitude, this 
means that power will be transferred from A 1 to A3 without 
ever going through A2. Note that the counter-intuitive order 
is maintained in the nonlinear optics case: first A2 and A3 are 
coupled, and the coupling between A 1 and A2 is introduced at 
a later point. 

In order to satisfy the adiabatic condition, it is required 
that the coupling between the desired eigenstate, \go), and the 
other eigenstates, |^±), is small compared to the difference 
between the effective wavenumbers of these states iQl 
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Figure 2. Normalized coupling coefficients of the two nonlinear pro- 
cesses with Gaussian modulation. Note that the 0)2 £03 coupling 
is (counter-intuitively) maximized before that of tOi o (B2. 



(7) 



This condition can be written in terms of 9 as 



d9_ 

Iz 



(8) 



where 9 is the adiabatically varying parameter, going from 
Oto k/2. 

Optimization of the atomic STIRAP process by temporally 
shaping (f), and also using time varying detunings (f), 
has been studied extensively fo - flal . Here we will use Gaus- 
sian modulation of the nonlinear coupling in space (along the 
interaction), which is convenient for analysis and illustrates 
the main features of STIRAP. For this purpose we replace the 
coupling coefficients in the matrix M of Eq. |2] with 



{z) 
K32 (z) 



Kne 



-(z-L/2-s)V»'' 
-{z-L/l+sf/w^ 



(9) 



where s and w are parameters which determine the locations 
of coupling maxima and the rate of the coupling variation, re- 
spectively. See fig. |2]for an illustration of these modulation 
functions. Note that the s parameter determines the coupling 
order: for < the intuitive order is obtained, while for s>Q 
the modulation is in the counter-intuitive order regime. Addi- 
tionally, as before, kij = icofkj/cojkA k*j^. 
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III. NUMERICAL SIMULATION UNDER IDEAL 
CONDITIONS 

A. Interaction in the Absence of Absorption 

In this section the main features of optical STIRAP will be 
demonstrated by numerical simulation, which solves Eq. [T] 
For this purpose, it will be assumed that both processes are 
simultaneously phase-matched, and that the desired modula- 
tion of the coupling coefficient kjj has been achieved. In this 
respect, the conditions are considered ideal. Section HV] will 
discuss a technique to imitate this situation using currently 
available technology. 

In the simulation, the two processes were SFG and 
DFG: (Ai = 3000«m) + (A,,, = SOOnw) (A2 = 631. 6nm) 
and (A2 = 631nm)- (A,„ lOOOnm) (A3 = 1714nm). The 
nonlinear medium was a 35mm long KTIOPOa (KTP) crys- 
tal with x^^^ («! , ; a>2) w (0)2, OJp, ; a>3) = 16.65/?»i/y 
i23]l . The coupling coefficients refractive indexes were cal- 
culated using the Sellmeier equation of Fradkin et. al iflTll 
and Emanueli et. al ifisll at a temperature of lOO^C. The 
modulation parameters were chosen to be s = 5mm and w = 
%mm. Note that for every i > the interaction was per- 
formed using the counter- intuitive order. The input intensity 
was IQOMW / cm^ and each of the two pumps had an intensity 
of 2GW /cm^. Figure [3] shows the resulting intensities of the 
interacting waves along the nonlinear crystal, with the inter- 
mediate wave intensity on a smaller scale in the inset. The 
input power is seen to be fully converted to the output wave, 
with some additional power from the first pump: if every pho- 
ton at Al is converted to a photon at A3, the intensity ratio is 
Ij/Ii = 0)3/0)1 = A1/A3. For the parameters used here, full 
conversion means 73 = 17 5MW/cm^. Furthermore, the inter- 
mediate power is at most 0.8% of the input power. These are 
exactly the characteristics of STIRAP. 

Interestingly, using the method of adiabatic elimination [01 
for the same set of parameters and under the same ideal con- 
ditions, the output intensity would reach only 73. 5MW /cm^. 
Clearly, the STIRAP analog method introduced here is more 
efficient, while maintaining the property of negligible inter- 
mediate frequency power 

For comparison, we repeated the same simulation with s — 
—5mm, i.e. with intuitive order The results are depicted in fig. 
m showing strong oscillations of all three waves intensities. 
Specifically, significant energy is present in the intermediate 
frequency, in contrast to the counter-intuitive result of fig. |3] 

Furthermore, we explored the conversion process dynamics 
as a function of the shift parameter s. In order to do this, the 
width parameter was kept constant atw — %mm while the shift 
parameter was varied from —20mm to +20mm. Efficient con- 
version from Al to A3 together with negligible (< 1%) power 
at A2 was obtained for 3.3mm <s < 6.6mm. In this shift inter- 
val, which is in the counter-intuitive regime, the overlap be- 
tween the modulation Gaussians of Eq. |9]is significant enough 
to facilitate STIRAP. It was also found that for values ofs < 0, 
i.e. in the intuitive regime, there is always high intermediate 
power somewhere along the crystal. 
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Figure 3. Numerical simulation of the intensities of the interacting 
waves along the nonlinear medium, with no intermediate wave ab- 
sorption, under ideal conditions and counter-intuitive order modula- 
tion {s = 5mm). The inset shows the intermediate wave intensity on 
a smaller scale. 
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Figure 4. Numerical simulation of the intensities of the interact- 
ing waves along the nonlinear medium, with no intermediate wave 
absorption, under ideal conditions and intuitive order modulation 
(s = —5 mm). 



B. Interaction in the Presence of Intermediate Frequency 
Absorption 

One of the key advantages of having a dark intermedi- 
ate state is that the crystal absorption in the intermediate 
wavelength does not change the conversion efficiency [Sj]. 
Here we consider the same interaction as the previous one, 
except that the intermediate frequency is in the ultraviolet 
absorption band of the KTP crystal. The two processes 
are (Ai = 2267«m) + {kp^ = 400nm) (A2 = 340nm) and 
(A2 = 340«m)- (Ap2 = ISOOnm) (A3 = 439.7«m). The ab- 
sorption coefficient at the intermediate wavelength is a2 — 
229. 9cm^' Ml9tl . Absorption is known to be accompanied by 
enhancement of the nonlinear susceptibility, and this is the 
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case here as well. However, this enhancement is not very high, 
since only one wavelength (A2) in each interaction is near an 
electronic transition resonance of the KTP crystal. We can 
calculate this enhancement using the density matrix quantum 
mechanical model of the linear and nonlinear susceptibilities 
lUt], which yields 
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where the primed co and 7 terms are resonant transition fre- 
quencies and full widths at half maximum, respectively, and 
COj and cOk are the two low frequencies out of the three in the 
TWM process. In the case under study, each of the two low 
frequencies in each of the two processes, namely Oi ,COpi, C0p2 
and 0)3, are far from any resonance. The last two terms of 
Eq. nn therefore do not contribute to the resonant enhance- 
ment of . The enhancement stems only from the first term, 
since it refers to the high frequency which is the intermediate 
frequency CO2 in each of the two processes considered here. 
From Eq. [TO]and[TT]is it therefore clear that in this case the 
second order susceptibility scales like the first order suscepti- 
bility. The enhancement factor for the wavelengths involved 
here can be calculated using the relation between the linear 
susceptibility and the refractive index, by making use of the 
Sellmeier equation of KTP: 



(1064nm+ 1064nm^- 532nm) 



A3) 



^(2) (1064nm+ 1064«m^ 532nw) 



n2 (A2) - 1 



(1064nm) (I064nm) - 1 



1.4 



(12) 



where we have used the 1064nm second harmonic generation 
value of as a reference. This enhancement is equivalent to 
increasing the intensity of each of the two pumps by a factor 
of 1.4^ « 2. Using these parameters together with the same 
input intensity, pump intensities and Gaussian modulation of 
the coupling coefficients as before, Eq. [T|was solved once 
again. The resulting intensities along the crystal are displayed 
in fig. |5] which shows that the conversion from Ai to A3 is still 
efficient (here full conversion means I3 = 515MW / crn^ due 
to contribution from the first pump) . The intermediate wave 
intensity, displayed in the inset, reaches at most 0.17% of the 
input intensity. These results reveal the true merit of avoiding 
significant power at the intermediate frequency, as it allows 
the conversion to remain efficient in spite of great absorption 
at A2. 



Figure 5. Numerical simulation of the intensities of the interacting 
waves along the nonlinear medium, with high intermediate wave ab- 
sorption and under ideal conditions. The inset shows the intermediate 
wave intensity on a smaller scale. 



IV. NUMERICAL SIMULATION USING 
PHASE-REVERSAL QUASI-PHASE-MATCHING 

In the preceding section it was assumed that both processes 
were phase matched and that the coupling coefficients were 
modulated as desired, without accounting for any mechanism 
to achieve these purposes. Now these assumptions will not be 
made. Instead, a technique called phase-reversal quasi-phase 
matching (PRQPM) [6] will be employed to provide phase- 
matching as well as to modulate the coupling coefficients. 
In PRQPM, binary functions are multiplied to produce a de- 
sired modulation, which has multiple Fourier components. 
This modulation can therefore provide QPM for multiple pro- 
cesses, i.e. processes with different phase-mismatches. For 
example, we can construct a product of two binary functions 
as follows: 



g (z) = sign 



sign 



-cos[7zDi)+cos[ —z 
—cos{7iD2) +COS ( 



(13) 



For each periodic term in the product, Dj is the duty cy- 
cle and Ay is the period. As always in QPM, we keep only 
the first order Fourier terms that contribute the desired phase- 
mismatches. If we choose Ai = A^i and A2 = Ak2, using sim- 
ple Fourier analysis yields 



gQPM (z) ~ — {2D2 - 1 ) sin {nD\ ) exp {±iAki ) + 
2 

- (2Di - l)sin(nD2)exp(±iAk2) (14) 
7t 

Evidently, this choice of the modulation periods provides 
QPM for the desired processes, while Di and D2 determine 
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Figure 7. Numerical simulation of the intensities of the interact- 
ing waves along the nonlinear medium, with no intermediate wave 
absorption and where the crystal is modulated using phase-reversal 
quasi-phase-matching. 



Figure 6. Normalized coupling coefficients of the two nonlinear pro- 
cesses with the modulation used in section ITVl The top panels show 
the PRQPM poling, which was used in the simulation of fig[7l at the 
beginning, center and end of the nonlinear crystal. The colored and 
white stripes represent domains with positive and negative nonlinear- 
ity. 



the the magnitude of the effective coupling coefficient for 
each process. Varying the two duty cycles along the crystal 
achieves the required modulation. 

The two cases regarded in section |III] with ideal phase- 
matching and coefficient modulation, are considered again in 
this section. In addition to using PRQPM for phase-matching 
and coefficient modulation, the calculation made here also 
takes into consideration technological restrictions on nonlin- 
ear modulation, and pumps depletion/amplification. We shall 
now assume that QPM is achieved by the wide-spread tech- 
nique of electric field poling of a ferroelectric crystal. 

Using QPM, especially when considering the technological 
limitations of this technique, makes it more difficult to satisfy 
the adiabaticity condition of Eq. |8] which essentially requires 
the changes in the coefficients to be very gradual. For this 
purpose, we first note that the spatial area where K32 is ris- 
ing while K'i2 ~ is going to waste, since = all along it 
(see fig. |2|l. In order to achieve adiabatic interaction, it is 
sufficient that at the start K"i2 ^ K32. We can therefore start 
the interaction with K32 (0) = max [K32 (z)], i.e. by choos- 
ing s — L/2, and use the additional space, which was previ- 
ously wasted, to make the coefficients gradients smaller, by 
increasing the modulation width w. For this reason, in the 
simulations conducted using PRQPM, the shift parameter was 
s = L/2 = 17.5mm and the width parameter was w = 25mm. 
The resulting normalized coefficients are displayed in fig. |6] 
along with an example of a corresponding PRQPM poling pat- 
tern at various locations. 

Furthermore, we note that Zukauskas et. al succeeded in 



fabricating domains as small as 345nm in KTP 120(1 . In the 
simulation conducted here, domains formed by Eq. [13] which 
were less than 350nm long, were concatenated with neighbor- 
ing domains until they reached this minimum length. In this 
manner modulation patterns that are attainable with current 
technology have been constructed. 

Simulation results for the first set of wavelengths, all of 
which are in the crystal's transparency spectral region, are dis- 
played in fig. |7] Comparing with the ideal case results shown 
in fig. [3] we see that the main features of the interaction have 
been preserved: high conversion from X\ to A3 is obtained, 
while the intermediate wavelength intensity remains relatively 
low. Still, there are two clear differences from the ideal case: 
the intermediate intensity reaches a higher value (9.15% of 
the input) and the intensities of all three waves fluctuate along 
the crystal. Both of these differences come from the use of 
PRQPM and the domain size limit: the oscillations are a well 
established property of QPM, coming from the non-phase- 
matched part of the interaction, i.e. higher Fourier orders of 
the modulation. The desired adiabatic variation is obtained on 
an average scale, which includes many domains. As a result, 
on a shorter scale, energy does get transferred to the interme- 
diate wave. The modulation periods being Ai = 18.2/xm and 
Al = 15.5/xm, the domain size limit makes it more difficult to 
satisfy the adiabatic condition, resulting in more energy trans- 
fer to A2. 

Simulation results for the second set of wavelengths, which 
include a highly absorbed intermediate wave, are depicted in 
fig. [8] While significant conversion from Ai to A3 took place 
in this case, it has some features which clearly distinguish it 
from the ideal case, in addition to the usual QPM oscillations 
noted above. First, the intermediate intensity reaches a value 
as high as 21.1% of the input intensity, where significant in- 
termediate power is only present in a confined spatial area. 
The reason for this phenomenon is that in this area the do- 
main size limitation causes loss of adiabaticity, resulting in 
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Figure 8. Numerical simulation of the intensities of the interact- 
ing waves along the nonlinear medium, with high intermediate wave 
absorption and where the crystal is modulated using phase-reversal 
quasi-phase-matching. 



energy transfer from A3 to A2, which is quickly absorbed by 
the crystal. Note that the modulation periods are Ai = 2.8/im 
and A2 = 2.1 jim, so small domains are abundant. Due to the 
loss of adiabaticity, the remainder of the interaction is also 
characterized by conversion from A3 to A2, albeit at a slower 
rate since little energy is present at the intermediate frequency 
at any point, due to the strong absorption. 

Optical STIRAP via PRQPM, in the presence of interme- 
diate frequency absorption, is ultimately limited by the ra- 
tio between the absorption length l^hs —Xj do. and the coher- 
ence lengths zf/'' = ;r/ ISkj , = 1 , 2 [ 1 ] . In order for any QPM 
scheme to have a significant effect, it is required that > Z^., 
otherwise an interacting wave will be absorbed before being 
affected by the modulation (for the case considered here, we 

have labs/lc^^ — 31.9 and labs/^?'^ = 41.1). Nevertheless, op- 
tical STIRAP via PRQPM still greatly reduces the conversion 
efficiency's sensitivity to intermediate frequency absorption, 
since considerably less energy is ever present in the interme- 
diate wave. Furthermore, no critical control of the modulation 
is required, as is the case for non-adiabatic QPM methods in 
the presence of absorption, which significantly limits the in- 
teraction length. 

An additional limitation of this modulation scheme is that 
in order to achieve independent control of each of the two pro- 
cesses via the duty cycles Di and D2, their phase-mismatches 
need to be sufficiently different from one another Otherwise, 
each of the two terms in Eq. [14] will have the same effect on 
both processes, and independent control would be lost. 



V. FUTURE PROSPECTS AND APPLICATIONS 

In the previous section limitations of PRQPM in face of 
current fabrication capabilities were discussed. Specifically, 
it was demonstrated that the minimum domain size limits 



this method's capability of satisfying the adiabaticity condi- 
tion when the dispersion is stronger, i.e. when the interacting 
wavelengths are relatively short (in the near UV). We con- 
sider it reasonable to assume that the current trend of ever 
decreasing available domain size will continue [20], making 
PRQPM feasible to processes involving short wavelengths, 
going into the UV. Note also that the domain size limitation 
can be somewhat circumvented, by relinquishing single pho- 
ton resonance A^i = Ak2 = and satisfying only two photon 
resonance A^i + Ak2 = 0. However, this will come at a cost of 
efficiency. 

Moreover, conversion efficiency can be further improved by 
employing techniques originally invented for atomic STIRAP 
113-16]. In these techniques the temporal dependence of the 
Rabi frequencies and detunings are controlled by pulse shap- 
ing and chirping mechanisms, respectively. Both of these con- 
trols have counter-parts in PRQPM optical STIRAP. As was 
shown here, the coupling coefficients can be controlled via the 
duty cycles of the PRQPM modulation function of Eq. [141 in 
analogy with pulse shaping of the Rabi frequencies. Addi- 
tionally, the effective phase-mismatch of each process can be 
modulated along the interaction by variations in the periods of 
the modulation (Ai and A2 in Eq. [T4l i. This is the counterpart 
of variable atomic detunings, obtained by chirping. 

An entirely different method for obtaining adiabatic in- 
teraction can involve anomalous dispersion phase-matching 
(ADPM) and tight beam focusing. In ADPM, absorption res- 
onances are deliberately brought forth, either by doping bulk 
nonlinear crystals with metals or ions 1I2TI1 . or by introducing 
chromophores in the fabrication process of polymer waveg- 
uides ]22]. These resonances influence the dispersion prop- 
erties of the nonlinear medium, thus enabling birefringent 
phase-matching of various processes. Nonlinear optics STI- 
RAP can be facilitated by using ADPM for phase-matching, 
while focusing the two pump beams at different locations 
along the crystal, achieving the modulation of the coupling 
coefficients simply by allowing diffraction to modulate pumps 
intensities (working with focused beams would invalidate the 
plane-wave approximation used in this work, however we be- 
lieve that this would not prevent our method from being ap- 
plied, and that it can be compensated by increasing pump 
power). Furthermore, the intermediate frequency can reside 
within the absorption band of the doped crystal without hin- 
dering the conversion efficiency. 

Two applications of nonlinear optics STIRAP will now be 
put forth. First, a straightforward application would simply 
be frequency conversion across large frequency ranges. By 
employing two SFG or two DFG processes, efficient conver- 
sion can be performed where the input and output frequencies 
are greatly separated from each other (e.g. mid-IR input and 
blue visible output, or vice versa). Second, we suggest that 
this method can be employed in the study of nonlinear crystal 
electronic structure and properties, by employing a SFG pro- 
cess followed by a DFG process. While Raman spectroscopy 
enables the study of molecular structure in the infrared ab- 
sorption band, experimental data on the electronic structure 
in the UV absorption band is scarce [|23ll24ll . Experiments to 
this affect rely mostly on nonlinear absorption and refraction 
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lHHH^, and are limited by linear absorption. UV electronic 
resonances will have a significant influence on nonlinear op- 
tics STIRAP processes, first via phase-matching and second 
via enhancement of the nonlinear coefficient. For example, 
by using a tunable pump or by tuning the crystal temperature, 
phase matching conditions can be controlled, thus allowing 
direct measurement of dispersion inside the UV absorption 
band (i.e. beyond the band gap energy). 

VI. CONCLUSION 

This work theoretically demonstrates a method for ef- 
ficient frequency conversion through an intermediate fre- 
quency, which never receives any significant amount of power 
throughout the interaction. This property was predicted an- 
alytically by analogy with atomic STIRAP, and confirmed 
with numerical simulations. It was shown that the absence 
of power at the intermediate frequency renders the conversion 



efficiency highly insensitive to absorption of the intermediate 
frequency, opening the way to conversion through absorptive 
bands in the UV. A technologically feasible method for carry- 
ing out such a process was proposed and numerically demon- 
strated to be effective. We suggest to apply this method in 
conversion between highly disparate frequencies, e.g. mid-IR 
to visible conversion. Another application is the study of the 
electronic structure and properties of nonlinear optical crys- 
tals, which is manifested in the UV absorption band, now 
made accessible in spite of linear absorption. 
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